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This paper presented two general quantum search algorithms. We derived the iterated 
formulas and the simpler approximate formulas and the precise formula for the amplitude 
in the desired state. A mathematical proof of Grover's algorithm being optimal among the 
algorithms with arbitrary phase rotations was given in this paper. This first reported the 
non-symmetric effects of different rotating angles, and gave the first-order approximate phase 
condition when rotating angles are different. 
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1 Introduction 



Shor reported his prime factoring algorithm in ||. Then Grover gave his quantum search 
algorithm In [§] Grover's quantum search algorithm consists of a sequence of unitary 

operations on a pure state. The algorithm is Q = —Iq WIt W, where li = I — 2\x)\(x\ and 
inverts the amplitude in the state \x), W is Walsh-Hadamard transformation. It would carry 
out repeated operations of Q, that is, ...{-I^WI^W^-I^WI^W^-I^WI^W)... 

= ...{-WI^ ) W)I^\-WI^ ) W)I^\-WI^ ) W)I^ ) where (-WI^W) is just the inversion- 
about-average operation. Therefore Grover's algorithm consists of alternating iteration of (—WI^W), 
1 The paper was supported by NSFC and partially by Lab. of computer science, ISCAS 
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inversion about the average, and J-f , inversion of the amplitude in the desired state |r). Then 
in J|| Grover showed how to replace Walsh-Hadamard transformation in his original algo- 
rithm with an arbitrary quantum mechanical operation, obtained the quantum search algorithm 
Q = -4 w) U- X I^U, where U is an arbitrary unitary operation and U 1 is the adjoint (the com- 
plex conjugate of the transpose) of U. Grover thinks that it leads to several new applications and 
broadens the scope for implementation. When the amplitudes are rotated by arbitrary phases, 
instead of being inverted, that is, Q = —I^U~ l I^U , where ii = I— (— e te + Long 
et al. first found 9 and (j) must satisfy a matching condition: 9 = (J), and derived an approximate 
formula for the amplitude in the desired state || and studied the effects of using arbitrary phases 
in amplitude amplification 0. In || the phase condition that tan((/?/2) = tan(^/2)(l — a) was 
presented. In [[H]] the recursion equations were used to study the quantum search algorithm. 

Here, this paper presented two general quantum search algorithms, derived the iterated 
formulas and the simpler approximate formulas of the amplitudes in the desired state. We 
showed that the amplitude in the desired state can be precisely written as a polynomial in 
(PX) for any quantum search algorithms which preserve a two-dimensional vector space. This 
is the first precise formula of the amplitude amplification in the desired state for the general 
quantum search algorithms including Grover' and Long et al.' and Hoyer's ones. The precise 
formula for the amplitude can help derive a precise phase condition. A mathematical proof was 
given in this paper for that Grover's algorithm is optimal among the quantum search algorithms 
with arbitrary phase rotations. We first found that the effects of rotating angles in the initial 
state and the desired state on the amplitude of the desired state are not symmetric. We also 
discussed the amplitude amplifications, and gave the first-order approximate phase condition 
when rotating angles are different. 

2 More general quantum search algorithm Q = —IjVI T U 

Let ijp = I — {—e ld + l)|sc)|(a;|. Grover studied the quantum search algorithm in [3]: Q = 
-I^U^I^U, where U is an arbitrary unitary operator and U 1 is the adjoint (the complex 
conjugate of the transpose) of U and Ix inverts the amplitude in the state \x). Generally let 
I x = I— ae lS \x )|( x\. Then I x is unitary if and only if (1— ae ld ){\— ae _,ff ) = 1. That is, a = 2 cos 9. 
Then I x = I — 2cos0e ie |x}IH/When e = Ix = jW_ If let j/ = j_ ( ae *e + i)| x )|( x | ) then 

I' x is unitary if and only if a = il.That is, l' x = Ix or I x 7T+9 \ Clearly I x w+B) = I x ^ I x and 
I x ^ = I - 2 cos %e ie \x)\(x\ and I x = I - 2cos9e i6 \x)\{x\ = I ( x +9) I { x ] = I x n \l x 8) ) 2 . Long et 
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al. studied the phase matching condition for the algorithm || Q = —I^U^I^U. 

Let's study the quantum search algorithm Q = —I^VI T U, where V and U are arbitrary 
unitary N x N matrices , where N = 2 n for n qubits, and I~ = I — 2 cos 9e lS \j )| ( 7I and 
I T =/-2cos0e^|r)|(r|. 

Since V and U are unitary, (VU) is also unitary, where (VU) is the product of V by U. 
Then by the definition (VU)^ 1 = [VU) + Jl]] which is the adjoint (the complex conjugate of 
the transpose) of (VU). Assume that (VU) is hermitian. By the definition (VU) = (VUy 
0, therefore (VU)(VU) = (VUVU) = I, where (VUVU) is the product of the matrices. In 
fact, given that (VU) is unitary then (VU) is hermitian if and only if (VUVU) = I. We give 
a proof as follows. Since (VU) is unitary and hermitian, (VU)(VU) = (VU)(VU) + = I. 
Conversely, from that (VUVU) = I, we obtain (VU)(VU) = (VU)(VU) + since VU is unitary, 
then (VU) = (VU) + . Therefore by the definition VU is hermitian. 

For Grover's algorithm V = U^ 1 clearly (VU) is hermitian. Here we give a sufficient con- 
dition in which (VU) is hermitian as follows. If V and U are hermitian and V and U are com- 
mutative, then (VU) is hermitian since by the definition V = V + and U = U + and (VU) + = 
u + v + = UV = VU. 

On the assumption that (VU) is hermitian we will show that Q preserves the four - di- 
mensional vector space spanned by I7), (VU)\j), V\t) and C/ _1 |r). After obtain V\t) and 
U \t), we apply V' 1 to V\ t) or U to U 1 |r), then obtain the desired state |r). Let's cal- 
culate the amplitude in the desired state V\r) and J7 -1 |t). We will use the notations in ^j, 
(t|C/| 7 ) = U TJ , (j\V\r) = V 1T , (rlV-ifr) = V* (t\(UV)\t) = (UV) TT , (j\(VU)\ 7 ) = (VU) 



After calculating Q\j) ,Q((VU)\j)) , Q(V\t)) and Q(U 1 |r)) we obtain the following express. 
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2.1 The iterated formula for the amplitude 

Next we derive the iterated formula for the amplitudes in the states | 7 ), (V?7)| 7 ),V|t) and 
[7 _1 |r) for Q k \^) after k operations of Q. 
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Let Q\ 1 ) = l 1 \ 1 )-{VU)\ 1 )+ Pl {V\r)), 
Q((VU)\ 7 )) = l 2 \ 1 )+p 2 (V\r), 
Q(V\r)) = hh)+P3(V\T))-U- 1 \r), 
Q(U~ 1 \t)) = U\ 1 )+ Pa (V\t)). 
Let Q k \j) = a k \j) + b k ((VU)\^)) + c k (V\r)) + d k (U~ l \r)), where a k ,b k , c k and d k are the 
amplitudes in the states I7), {{VU)\^)), (V|r)) and (U~ x \t)), respectively. Then Q fc+1 |7) = 
Q(Q k \ 1 )) = a k Q\i) + b k Q((VU)\j)) + c k Q(V\r)) + 4Q(C/- 1 |r))- Th en we obtain the following 
iterated formula:a fc+1 = ha k + hh + hc k + hd k , b k+1 = -a k , c k+1 = p\_a k + p 2 b k +p?,c k +p^d k 

> <4+l = —Ck ■ 

2.2 The first-order approximate formula for the amplitude 

From the iterated formula we will approximate c k+ \ using the first-order Taylor formula of 
U T -y. That is, in c k+ \ we only keep the first-order of t/ T7 and omit the high order of U T1 . Let 
Vy T = U*, (VU)yy = (UV) TT = to approximate c k+ \. Let's give a brief justification as 
follows. Clearly if V = ±U~ X , then (VU) = ±1, then {VU) yy = {UV) TT = ±1. Therefore 
assume that V / ±U~ 1 . Let Wij be any term of a unitary matrix W. Then by the definition 
WW + = W + W = 10, J2 w kiW k j = J2 w ikWjk = %> where §ij is the Kronecker deta and w k j 

k k 

is the complex conjugate of w k j. That is, each row(column) of W has length one and any two 
rows(columns) of W are orthogonal. Clearly \wij\ < 1. If \wij\ = 1, then all terms of the ith row 
and the jth column of W are zero except the ij term. The results above hold for V, U, (VU) 
and (UV) since they are unitary. Specially |(V£/) 77 | < 1 and |(C/V) Tr | < 1. 

If W is assumed to be hermitian |Q| 3 then by the definition W = W + ,tlia.t is, w^j = Wji, then 
all the diagonal terms ii of W are real. Since V is unitary and (VU) is assumed to be hermitian, 
it is easy to show that (UV) also is hermitian. By the definition (VU) = (VU) + = U + V + , then 
V + (VU)V = V + (U + V + )V, UV = V + U + = (UV)+, therefore (UV) also is hermitian. So all 
the diagonal terms of (VU) and (UV) are real. Specially (V?7) 77 and (UV) TT are real. 

Next we will give a sufficient condition in which all the diagonal terms of (VU) and (UV) 
are zero. For the detail please see the appendix 4. 

If V is unitary and each Vij of the block form of V is of the form ^ ^ ^ , and U = V + P. 

Then U is unitary, (VU) and (UV) are unitary and hermitian, and all the diagonal terms of 
(VU) and (UV) are zero. For the definition of P please see the appendix 4. 

Now let us make the first-order approximate formula for the amplitude. Since p\ and p 2 
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contain the factor U T -y, we have to approximate with the expression without t/ r7 .From Q\j) 

0, c 2 = -2 cos (pe^Urj, 
0, k is odd; 

(-l) m e i2me , k = 2m. 



and It) °i = 0, fei = —1, c\ = 2 cos ^e l ^J7 T7 , (ii = 0; 02 = 
c?2 = — 2 cos (j)e l ^U ri . By induction we obtain a^+i = — e l2e afc_i, 



c k+ i = 2cos(/>e^(a fc -a fc _i)f/ r7 +(l-2cos0e^)c fc _i, c fc+ i = < 



(-l) m 2cos^[/ T7 Y. <? (m ~ l) 5\ 



1=0 



where a = 2 cos Oe — 1 



J20 



Case l,cr = 5.That is, 6 - 
terminology. 

-l) m (fe + 2)cos0e^ m C/ r7 , 
-l) m+1 (k + 1) cos (\>e^a m U, 



(_!)m+i 2 cos «^*7 T7 £ (j^-')^, 
z=o 

(5 = 2 cos - 1 = e i2?i . 

, which is called phase matching condition by Long et al.'s 



In the case Ck+\ = 
l c fe+i| " 



T~fl 
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k 



2m; 

2m + 1, 



(k + 2)| cos </>||£/ T7 |, k = 2m; 
(k + l)|cos0||C/ T7 |, k = 2m + \. 
Therefore for arbitrary unitary operators V and J7 , when (VU) is hermitian, Q = —I^VI T U 

can be used to construct a quantum search algorithm that succeeds with certainty except that 

I 7 = I or I T = /. The number of iterations is almost 1/ |cos <^>| |J7 T7 | to reach the desired state 

|r) from the initial state I7). 



Specially when a = 5 = l,that is, 9 
becomes — I^VI^U, and 
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We can derive the above independently, please see the appendix 2. 

|„, I _i. 9 |„™ j.1 irr I sin(m+l)(0-</>) 
|C fe+ i| - ^ |COS0| \U T1 \ sin ( e _^) — 

(& + 2)|[/ T7 | j cos 4>\, k = 2m; 



Case 2, <r / 5, that is, / 
2m + 1. Clearly lim |cfc + i| = 



(k + l)|C/r 7 || cos 0|, fc = 2m + l. 



, where = 2m or 
It means when almost \9 — <j>\ < 



2\ cos 0||f7 r7 |, Q can be used as a quantum search algorithm that succeeds with certainty though 
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3 The precise formula for the amplitude for the general algo- 
rithms which preserve a two-dimensional vector space and the 
non-symmetric effects of different rotating angles 



Though in the algorithm Q = —I^VI T U in the section 1 above let V = U then obtain the 
algorithm in this section Q = —I^U~ l I T U , note that in the section 1 when we derived the 
iterated and approximate formulas of the amplitudes in the state V\t) we assumed V 1T = f/* , 
(VU)yy = (UV) TT = 0. And in the section 1 Q preserves the four - dimensional vector space 
spanned by (7), (VU)\j), V\t) and ?7 _1 |r).In this section we will show that Q preserves the 
two-dimensional vector space spanned by (7) and £/ _1 |t), clearly V JT , (FJ7) 77 and (UV) TT don't 
appear. So the results in this section can not be obtained from the ones in the section 1 by 
simply letting V = U~ l . So for the algorithm in this section it is necessary to derive its iterated 
formula, precise one and approximate one of the amplitude in the desired state and discuss its 
amplitude amplification. 

Let's study the algorithm Q = —IjU^ItU, where 7 7 = I — 2 cos 9e ie \^j )|( 7) and I T = 
I — 2cos0e^|r)|(r|. When = <j> = 0, it reduces to Grover's algorithm. 

After calculating, we obtain the following expression. For detailed derivation, please see the 
appendix 3. 



where a = -(1 - 2 cos 0e id + 4cos 6e ie cos 0e^|£/ r7 | 2 ), (3 = 2 cos 0e^f7 r7 , A = 2cos0e^(l - 
2 cos 4>e lr ^)U*y, 5 = 2coscf)e 1 ^ — 1. Clearly the present algorithm Q like Grover's and Long et 
al.'s algorithms preserves the vector space spanned by I7) and ?7 _1 |t). 
For Grover's algorithm a = 1 - 4\U T1 \ 2 ,f3 = 2C/ r7 , A = -2U* 5 = 1. 



For Long et al.'s algorithm a = -e ie - ( -e id + l)(-e^ + l)|^ r7 | 2 , P = {-e i(p + l)U TJ , A = ( 
_ e i0 + i) e i<t>u* 5 = -e**. 



Hoyer's algorithm preserves a two-dimensional vector space spanned by another two states, 



where a = -{(1 - e**)a + e**},p = (1 - e^^E^/T^ae^ , A = (1 - e^^E^/T^a, S = {(1 - 
e**)a - l}e^. 



The process of deriving the following iterated formula and precise formula for the amplitude 
is the same for the present algorithm and Grover' and Long et al.'s and Hoyer's algorithms and 
any other quantum search algorithms which preserve a two-dimensional vector space, this is 
because the derivation is not concerned in the contents of a, (3, A and S. 
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3.1 The iterated formula for the amplitude 

Let's derive the iterated formula for the amplitude in the desired state [7 _1 |r) of Q k \^)- Let 

Q| 7 > = «|7> + PiU-'lr}), Q{U- l \r)) = A| 7 ) + ^V))- Then 

Q 2 | 7 > = (a 2 + /?A)| 7 ) +/?(« + ^(CZ-V)), (1) 

Q3| 7 ) = (a 3 + £A(2a + 5))\j) + /3(a 2 + a5 + 5 2 + /3A)(?7- 1 |t))v(2) 

Q 4 |7> = (a 4 + /3A(3a 2 + 2a5 + 5 2 ) + (/?A) 2 )| 7 > 

+/3(a 3 + a 2 5 + a5 2 + 6 3 + 2(a + <5)/3A)([/- 1 |r))...(3) 

Let Q fc | 7 ) = a/c| 7 ) + &fc(t/" — 1 |t)), where a& and b k are the amplitudes in the initial state | 7 ) 
and the desired state U~ x \t), respectively. Then Q fc+1 | 7 ) = Q{Q k \l)) = a k Q\j) + b k Q(U~ x \T)) 

= (aa k + A6 fc )| 7 ) + (f3a k + <56 fc ) 1 t> ) ,clearly a k+ ± = (aa k + Xb k ) and b k+ i = ([3a k + 5b k ). 

It is also the iterated formula for amplitude amplification for Grover's and Long et al.'s and 
Hoyer's algorithms and any other quantum search algorithm which preserves a two-dimensional 
vector space. Clearly it does not need to computer Q k {U~ l \r)) to derive the iterated formula. 

3.2 The precise formula for the amplitude 

From the iterated formula above by induction a k and b k can be precisely written as the following 
polynomial in (/3A), respectively. Let [x] be the greatest integer which is or less than x. 
b k = P(c k0 + c fcl (/3A) + c k2 (f3X) 2 + ... + c k[{k _ 1)/2] (PX)^- 1 ^), (4), 

where c kj = ^^x-^^ 1 ^^ = ^ ^ ) ( ^f 1 ). 

a k = a k + d kl {(3\) + d k2 (PX) 2 + ... + d k[k/2] (pX)^ (5), 

where d kj = J£ J k %_ 2j _ n) a n 5 k - 2 ^ , and ig> = ( ^ ) ( ^ ) .Note that ( * J = 

l,for any n > 0. 

For Grover's algorithm (/3A) = a — 1, J = 1. Then b k can be written as b k = (3r k , where r k 
is real. 

Here it does not need to computer Q k (U~ 1 \t)) to derive the precise formula above. Clearly 
it will save much more time to computer amplitude amplification in the desired state using the 
precise formula for b k than the kth power of the matrix which represents the operator Q. The 
formula is also the precise formula for the amplitude amplification for Grover's and Long et 
al.'s and Hoyer's algorithms and any other quantum search algorithm which preserves a two- 
dimensional vector space. From the precise formula for b k , it is not hard to see the phase 
conditions in || are only sufficient. 
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For the detailed derivation of bk, please see appendix 3. 

3.3 The non-symmetric effects of different rotating angles 

Let us study the non-symmetric effects of different rotating angles of the initial state (7) and 
the desired state 1 |r) on the amplitude in the desired state. Clearly the norm of bk contains 
the factor \@\ which is 2| cos^| |f/ r7 |.It is not hard to see that the effect of 4> on the amplitude 
in the desired state is greater than that of 9 when 9 / 4>.It means that when 9 / <p the 
effects of 9 and <p on the amplitude in the desired state are not symmetric. For example, when 
4> = vr/2, then I T = I,Q = —I^U^IrU = —I-yU~ l U = —Ij,bk = 0, it means that Q does 
nothing in the amplitude in the desired state except that it rotates the phase of the initial 
state | 7 ) by angle vr + 29. When 9 = it/2, then J 7 = I,Q = -I^U^UU = -XJ- x l T U ,Q k \i) = 

k—l 

(-l) fc | 7 ) + (/3 £ (-l) fe - 1 - i (5 i )C/" 1 |r)> \h\ = 2|£/ r7 ||sinfc(7r/2-0)|. The result in this subsection 

i=0 

is also true for Long et al.'s algorithm. 

3.4 The first-order approximate formula for the amplitude 

Long et al. used many transformations and approximate operations to derive the approximate 
formula for the amplitude amplification [4]. Next let's derive only by induction the first-order 
approximate formula for the amplitude bk in the desired state U~ 1 \t) using the iterated formula 
above though it is easy to derive it from the polynomial in ((3X) of a,k and bk, please see (4) and 
(5)above. We only keep the first order of U T1 in the amplitude in the state C/ _1 |r). In bk+\, (3 
contains the factor U T1 , so ak should be approximated with the express without U T1 ; 8 does not 
contain the term U T1 , so 5bk only contains the first order of U T1 provided that bk only contains 
the first order of U T1 . In afc +1 , since A contains the factor U* and bk contains the factor U T1 , 
Xbk must contain |C/ T7 | 2 and is omitted. Therefore a^ + i should be approximated by a ak, that 
is, dk+i = a ak, to make a^+i not contain factor U TJ . Let's see how to approximate Q k \"f) by 
induction. Clearly 

Q2 h) = a 2 h)+(3{(J + g){u -l lT) y 

k—l 

Assume that = a k \j) + 0{ £ <7 fc - 1 - < <J i )(E/- 1 |r)). Then Q k+1 \-f) = a k+1 \j) + ((3a k + 

i=0 

k—l k k—l 

5(3 £ a k - 1 - i 8 i ){U- 1 \T)) = (J k+1 h)+P(E a^S^iU- 1 ^)). Then by induction b k =(3 £ a k - 1 ~ i 8 i . 

i=0 i=0 i=0 
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It is easy to verify that the approximate formula is just the first-order of U T1 in the precise for- 
mula for bk- Please see the appendix 3. Clearly it does not need to computer Q k (?7 _1 |t)) to 
derive the approximate formula above. 

3.5 A mathematical proof of Grover's algorithm being optimal among the 
algorithms with arbitrary phase rotations 

Next let's use the approximate formula to study its amplitude amplification and prove that 

fc-l 

Grover's algorithm is optimal. Clearly \[3 & k ~ l ~ l li l \ < 2k \U T j\.We will prove that for any 

i=0 

case the norm of amplitude in desired state U~ 1 \t) after k operations of Q is less than 2k \U T y\ 
except Grover's algorithm. 

Case 1, a = S, that is, 9 = cj). In the case the amplitude b k = 2k cos (^e % ^a k ~ x U T ^, \b^\ = 
2k\ cos</>| |?7 r7 |. When cos(/> 7^ 0,that is, i~ = I T 7^ I, we obtain a quantum search algorithm that 
succeeds with certainty. 

Case 1.1. When a = 5 = 1, that is, 9 = 4> = 0,that is Grover's algorithm, 6& = 

fc-i 

/3 cr k ~ 1 ~ t 5 l = 2kU T1 , (please see the appendix 1 to check it), the norm of the amplitude 

i=0 

\bk\ = 2k \U T y\. Let 2k \U T1 \ = 1, we obtained the maximum probability for the desired state 
U~ 1 \t), in the case the optimal number of operations of Q is 1/2|C/ T7 | , when \U Try \ is taken as 
l/y/~N, the number is yN /2. Here the optimal number is less than Grover's optimal number 
7r/4|f7 r7 | evaluated in Q. The optimal number of iteration steps obtained by Long et al. || is 
7r/4/3,where \U T1 \ = sin/3,so the optimal number also is almost 7r/4|J7 T7 |. Please see the \bk\ in 
the table 1. 

Table 1. _ 

N VN/2 U t1 = 1/y/N k \b k \ 

100 5 0.1 6 0.9375 

400 10 0.05 12 0.9334 

625 12 0.04 14 0.9010 

900 15 1/30 17 0.9064 

k— 1 

Case 1.2. When a = S + l,clearly \(3 £ <r k - 1 - i S i \ < 2k \U T7 \. 

i=0 

Case 2, a ^ 5, that is, 9 ^ 0. Clearly £ a* -1- *** = ^E? ■ \h\ = \P E cr*- 1 "^ = 

i=o ' ' i=0 



2| cos ^H^ly \™ls2(Lt) = 2 I cos <i>\\ U T-f\\ s 7m{e-<p) I • In the case a ^ 008 - 4>) ± l M 

the induction on fc, it is not hard to prove that |~^j^r~^l < > 1)j therefore = 
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ft— 1 

\(3 £ v*' 1 "^^ < 2k \U T1 \. 

i=0 



From the cases 1 and 2, for Grover's algorithm in the first-order approximate \b k \ = 2k \U T ^\; 

for other cases \b k \ < 2k \U T ~\. It proved that Grover's algorithm is an optimal one with the 

f 2k\ cos (f)\\U T ^ a ~ A - 



form —I^U 1 I T U. Fvom above \b k \ = 
formula is simpler than Long et al.'s one B 



2 | cos 4>\ \U T1 \ 



sin fc(#— </>) 



. clearly the approximate 



3.6 The first-order approximate phase condition with different rotating an- 
gles 

In H Long et al. studied the effects of imperfect phase inversion. In || Peter Hoyer thinks if 

9 ^ (f> and \ip — 4>\ < c/y/N for some approximate constant c then the marked state can still 

be found by Long et al.'s algorithm with high probability. Here we will give the first-order 

approximate phase condition that \9 — (j>\ < 2\ cos ^>||f/ T7 | and deduce when 9 ^ 4> an d and (p 

satisfy the condition then the desired state can still be found by the present algorithm. 

Let Max k \bk\ be the maximal \b k \ for any k. When ^ <f> and \9 — cf>\ is small, Max k \b k \ = 

2\ cos \U Try \/\9 — It means that Max k \b k \ is the inverse ratio of \9 — <j)\ when \9 — <j)\ is 

small. When cos0 / 0,let \9 - <j>\ = 2l\ cos0||C7 r7 |. Then Max k \b k \ = l/l. When < / < l,that 

is, |6> — 0| < 2| cos 0||f7 T7 |, almost Max k \bk\ = 1; when 1 < Z,that is, \9 — <fi\ > 2| cos (/>||[/ T7 |, 

Max k \b k \ < l/l < 1. Please see the following table 2. In the table 2 let U rl = 1/y/N and (j) = 0, 

where N = 100, C/ T7 = 0.1. The experiments were done on IBM PC using MATLAB. 

The table 2. 
8 k \b k \ 

0.01 7 0.9899 

0.02 8 0.9994 

0.03 8 0.9930 

0.04 100 0.9861 

0.05 100 0.9525 

When a ^ 5, clearly lim 2| cos (f>\ \ U T1 \ \ m ^^gZ^ I = 2k\ cos 0||f7 T7 |. When = 0, the limitation 
is 2k \U T ~\, which is just Grover's algorithm. 
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For Grover's algorithm, Q = — 7 7 £/ 1 Ir U, 
Q\i) = (1 - 4|*7 T7 | 2 )| 7 > + 2U t1 U~ 1 \t) 

Q 2 It> = Q{Q\i)) = ((i - 4|f/ T7 | 2 ) 2 - 4|c/ T7 | 2 )| 7 ) + (4C/ r7 - 8[/ T7 |£/ T7 | 2 )(£/- 1 |t»- 
Q 3 | 7 ) = Q(Q 2 | 7 )) 

= ((1 - 4|f7 T7 | 2 ) 3 - 12|f/ r7 | 2 + 32|[/ T7 | 4 )| 7 ) + (6C/ T7 - 32[/ T7 |£/ T7 | 2 + 32?7 T7 |[/ T7 | 4 )(^- 1 | T )) 



Appendix 2 

Q = -I^VI^U. And note that ( 7 | 7 ) = 1, (r|C/| 7 ) = C/ T7 , ( 7 |C/ _1 |t) = £^ 7 ( Since 
U- 1 = U* and ( 7 |C/*|r) = (r|C/| 7 )* = 17^). Then 

Q| 7 ) = -/ 7 y/ T ?7| 7 ) = -(/-2| 7 )( 7 |)y(/-2|r)(r|)f/| 7 ) 

= -(J-2| 7 )( 7 |)FC/| 7 )+2(J-2| 7 )( 7 |)F|r)(r|C/| 7 ) 

= -(^)| 7 ) + 2| 7 )( 7 |yC/| 7 ) + 2?7 T7 (y|T))-4C/ T7 | 7 )( 7 |y|r) 

= -(W7)| 7 > + (2(VU)^ - 4*7 T7 VV)| 7 > + 2U TJ (V\r)). 

Q(U- 1 \t)) = -LfVIrUiU- 1 ^)) = -I 7 VI t \t) = I 7 V\t) 

= (I-2\ 1 )( 1 \)V\t) = V\t)-2V 1t \ 1 ). 

Q(VU\ 1 )) = -I 7 VI t U(VU\ 1 )) = -(I-2\ 1 )(i\)V(I-2\t)(t\)U(VU\ 1 )) 
= -(/ - 2| 7 )< 7 |)((y[/W7)| 7 > - 2V\t)(t\UVU\ 7 )) 

= -(VUVU)\j) + 2V\t)(t\UVU\j) + 2\j)(j\VUVU\j) - 4\j)(j\V\t)\(t\UVU\j) 

When VU is hermitian, then VUVU = 7, VUV = U' 1 and UVU = V~ l . 

Then Q(VU\ 7 }) = -| 7 ) + 2V; t (V\t)) + 2| 7 ) - 4V yr V* T h) = (1 - 4|y 7T | 2 )| 7 ) + 2V; t (V\t)). 

Q{V\t)) = -I,VI t U(V\t)) = -(I - 2\ 7 )( 7 \)V(I - 2\t)(t\)U(V\t)) 

= -{I -2\ 1 ){ 1 \){{VUV)\t) -2V\t){t\UV\t)) 

= -{(VUV)\t) - 2V\t)(t\UV\t) - 2\ 1 )( 1 \VUV\t) + 4| 7 )( 7 |V r |r)(r|?7V r |r)} 
Then Q(V\t)) = -U'^r) + 2{UV) tt (V\t)) + {2U* T1 - 4V^(UV) TT )\ 7 ) 
Appendix 3 

The algorithm Q = -I^U^IrU, where 7 7 = 7-2cos6»e* e | 7 )|( 7 | and I T = 7-2cos^e i *|r)|(r| 
Q| 7 ) = -IJJ-H T U\i) = -7 7 C/- 1 (7 - 2 cos <^|t ) | ( r\)U\i) 
= -7 7 (C/- 1 f/| 7 ) -2cos<pe i4 'U- 1 \T}\(T\U\"f}) 
= -7 7 (| 7 ) - 2cos<^C/ T7 (C/- 1 |t))) 

= -((7 - 2 cos 6e i9 \-f )|( 7 |)| 7 ) - 2 cos 4>e^U T1 (I - 2 cos fle i9 | 7 )| ( 7 |)(*7 _1 |t))) 
= -(1 -2cos0e i9 + 4cos^ e cos0e^|[/ T7 | 2 )| 7 ) + 2cos0e^?7 T7 (^ 1 |T))- 
Q([/ _1 |r)) = -7 7 c/" 1 7 t ?7([/- 1 |t)) = -7 7 ?7" 1 7 t |t) 
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= -J 7 tf _1 (J- 2cos0e^|r)|(r|)|r) 

= -I^U-^t) - 2cos<j)e i<t> U~ 1 \T)\{T\T)) 

= -J 7 (l - 2cos0e^)(?7- 1 |T)) 

= -(1 - 2cos0e^)(/ - 2cos#e ie |7)|(7|)([/- 1 |T» 

= -(l-2cos0e^)([/" 1 |r) - 2cos0e if? |7)|(7|[/~ 1 |T» 

= -(1 -2cos0e^)(C/^ 1 |r) - 2 cos Oe^U*^)) 

= (2cos0e^ - l)(^ 1 |r)) + 2cos6>e* e (l - 2cos<£e^)[/"* 7 |7) 

The following is the detailed derivation of the precise formula for the amplitude amplification. 
Let Q k \"f) = dk\l) + frfc(k r ~ 1 | r ))> where and b k are the amplitudes in the state I7) and the 
desired state ?7 _1 |r), respectively. Then 

6 5 = /3(( a 4 + a 3 5 + a 2 5 2 + a5 3 + 5 4 ) + (3a 2 + 4a<5 + 35 2 )/3A + (/3A) 2 ) 
a 5 = a 5 + (4a 3 + 3a 2 <5 + 2a5 2 + <5 3 )/?A + (3a + 2<5)(/3A) 2 

b 6 = /3((a 5 + a 4 5 + a 3 <5 2 + a 2 <5 3 + a5 4 + J 5 ) +(4a 3 + 6a 2 <5 + 6a5 2 + 4<5 3 )/3A + 3(a + <5)(/3A) 2 ) 
a 6 = a 6 + (5a 4 + 4a 3 <J + 3a 2 5 2 + 2a5 3 + 6 4 )(3X +(6a 2 + 6a<5 + 3<5 2 )(/?A) 2 + (/3A) 3 

From the iterated formula it is not hard to by induction show that a k and b k can be written 
as the following polynomial in (/3A), respectively. 

b k = (3(c k0 + c kl (f3\) + c k2 {f3\) 2 + ... + c k[{k _ l)/2] {[5\f k -^W), , 

where c kj = £ lgL-:y^)° B **~ 1 ~ !y ~ B - 

n=k-l-2j y J ' 

a k = a k + d kl (J3\) + d k2 (p\) 2 + ... + 4[fc/2](/?A) [fc/2] , 

where 4,= E t[%_ 2j _ n) a n 5 k - 2 ^ . 

n=k-2j y J ' 

l 3 ki in the coefficients of (5\ in 63, 64, 65, b$, 67 constitute the following pyramid. 
The table 3. 

1 h 
2 2 64 
3 4 3 65 
4 6 6 4 fe 6 
5 8 9 8 5 b 7 
Note that the diagonal elements 1,2,3,4,5,..., can be represented by ({), ( 2 ), ( 3 ), ( 4 ), (f)...; the 
diagonal elements 1,2,3,4,... times 2 are 2,4,6,8,..., , respectively,...; then they times 3 are 3,6,9,..., 
respectively;... 
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For example, 4o\ 4V' 4^> 4^ and 4^ i* 1 the coefficient C71 of /3A in 67 can be represented 
by (})(!), (?)(!), (?)(?), (f)(?) and (f)(|) respectively. 

l° ki in the coefficients of (/3A) 2 in 65, be, 67, 6s constitute the following pyramid. 
The table 4. 

1 65 
3 3 6 6 

6 9 6 b 7 

10 18 18 10 bg 
Note that the diagonal elements 1,3,6,10,..., can be represented by (2), (2), (2), (Dv > then 
the diagonal elements 1,3,6,... times 3 are 3,9,18,..., respectively; then they times 6 are 6,18,..., 
respectively;... 

For example, 4o\ 4? an< ^ 4f* ^ n ^ e coefficient C72 of ((3X) 2 in 67 can be also represented by 
(!)(!), and respectively. 

l 3 ki in the coefficients of (/3A) 3 in 67, 6g, 69, 610 constitute the following pyramid. 
The table 5. 

1 67 
4 4 6 8 

10 16 10 69 

20 40 40 20 610 

Note that the diagonal elements 1,4,10,20,... can be represented by(jj), (3), (3), (3),...; the 
diagonal elements 1,4,10,... times 4 are 4,16,40,..., respectively; then they times 10 are 10, 40,..., 
respectively;... 

For example, 1\q , l\ 01 , l 102 and l 103 in the coefficient C103 of (/3A) 3 in b\o can be also repre- 
sented by (|)(|), and (|)(|) respectively. 

Generally l 3 ki in the coefficients of {(3X) k in 62*1+1, &2fc+2, &2fc+3i &2fc+4, &2fc+5 constitute the 
following pyramid. 

The table 6. 

& k + 3 ) (t^C 2 ) ( k k + 2 )(1 +1 ) ( k k + 3 )C k ) "2 fc + 4 

OC 4 ) ( k k +1 )C 3 ) C 2 )C 2 ) ( k k +3 )( k k +1 ) C 4 )0 6 2 fc+ 5 

Therefore we can conclude and prove by induction that l^) = ^ ^ ^ ^ ^ V ^ ^ .So we 

obtain the precise formula for the amplitude b^ in the desired state after operations of the 
algorithm Q. 
Appendix 4 



13 



The sufficient condition in which (VU) and (UV) are hermitian and all the diagonal terms 
of (VU) and (UV) are zero 

Let's define the NxN matrix P. Let pij be any terms of P, where P(2i-i)(2j) = 1) P(2i)(2i-i) = 
1 and pij = otherwise. Clearly P is unitary and hermitian, p 2 = 1, and all the diagonal terms 
of P are zero. It is not hard to see that P has the block form which is diag{Pi, P2, ...,P m }, 

where m = N/2 and Pi is ^ ^ J ^ for i = 1, 2, m. 

It is easy to verify that (PV) means to interchange lines (2k — 1) and (2k) of V and (VP) 
means to interchange columns (2k — 1) and (2k) of V, where k = 1,2, ...,N/2. 

Given that V is unitary. Then if £/ = F+P then U is unitary and (VC) = V(V + P) = 
(VV+)P = P. AS well if V = U + P then (UV) = (UU + )P = P. From U = V+P, obtain 
U + = PV. From F = U + P, obtain V = PVP. Then we conclude the following sufficient 
condition. 

Sufficient condition (Version 1). If V is unitary, V = PVP and U = V + P, then (VU) and 
(UV) are unitary and hermitian, and all the diagonal terms of (VU) and (UV) are zero. 

Proof. From the discussion above (VU) = P. Since V = PVP, (UV) = (V + P)(PVP) = 
y+p 2 yp = (y + V)P = P. By the property of P, clearly the lemma holds. 

/ V u V 12 .. V14 \ 
V21 V22 •• V24 



Let V has the block form 



submatrix. 



, where m = N/2 and each Vij is 2x2 



V V m i V m 2 .. V mn J 



It is not hard to see that V = PVP if and only if Vij is of the form ( ° ] for 1 < i, j < m. 



b 



>ij 

From this we conclude the version 2 of the sufficient condition. 
The sufficient condition (Version 2) 

If V is unitary and each Vij of the block form of V is of the form ^ ^ ^ , and U = V + P. 

Then (VU) and (i/V) are unitary and hermitian, and all the diagonal terms of (VU) and (UV) 
are zero. 
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